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We derive the time-independent Gross-Pitaevskii equation at zero temperature for condensed 
bosons, which form as bound-fermion pairs when the mutual fermionic attractive interaction is 
sufficiently strong, from the strong-coupling limit of the Bogoliubov-de Gennes equations that de- 
scribe superfiuid fermions in the presence of an external potential. Three-body corrections to the 
Gross-Pitaevskii equation are also obtained by our approach. Our results are relevant to the recent 
advances with ultracold fermionic atoms in a trap. 
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Evolution from superfiuid fermions to condensed com- 
posite bosons appears on the verge of being experimen- 
tally realized in terms of dilute ultracold fermionic atoms 
in a trap, although several experimental difficulties still 
remain to be overcome 1 . Dilute condensed bosons in 
harmonic traps, in particular, have already been stud- 
ied extensively 2 . From a theoretical point of view, their 
macroscopic properties (such as the density profile) have 
been described by the time-independent Gross-Pitaevskii 
(GP) equation for the condensate wave function 3 (see 
also Refs. 4 and 5). Superfiuid fermions in the presence 
of a spatially varying external potential, on the other 
hand, are usually described in terms of the Bogoliubov-de 
Gennes (BdG) equations 6 , where a spatially varying gap 
function is obtained from a set of two-component fermion 
wave functions. Several problems (such as the calcula- 
tion of the Josephson current ) have been approached in 
terms of the BdG equations for superfiuid fermions. 

Given the possible experimental connection between 
the properties of superfiuid fermions with a strong mu- 
tual attraction and of condensed bosons, a correspond- 
ing theoretical connection between these two (GP and 
BdG) approaches seems appropriate at this time. In this 
way, one could even explore the interesting intermediate- 
coupling (crossover) region, where neither the fermionic 
nor the bosonic properties are fully realized. Such 
a connection might be also relevant for the study of 
high-temperature (cuprate) superconductors, for which 
the coupling strength of the superconducting interaction 
seems to be intermediate between weak coupling (with 
largely overlapping Cooper pairs) and strong coupling 
(with nonoverlapping composite bosons) . 

In this paper, we establish this connection by showing 
that the time-independent GP equation at zero temper- 
ature can be obtained as the strong-coupling limit (to 
be specified below) of the BdG equations. This result 
thus shows that the BdG equations, originally conceived 
for weak coupling, are also appropriate to describe the 
strong-coupling limit, whereby the fermionic gap func- 



tion is suitably mapped onto the condensate wave func- 
tion of the GP equation. Our derivation resembles closely 
Gorkov's microscopic derivation of the Ginzburg-Landau 
equations from BCS theory near the superconducting 
critical temperature in the weak-coupling limit 8 . Our 
approach further enables us to obtain high-order correc- 
tions to the GP equation; in particular, the three-body 
correction is here explicitly derived. 

We begin by considering the BdG equations for su- 
perfiuid fermions in the presence of a spatially varying 
external potential V(r): 
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Here 

H{r) = ~£ + v{r) * (2) 

is the single-particle Hamiltonian reckoned on the 
(fermionic) chemical potential fi (m being the fermionic 
mass and h = 1 throughout), while 

A(r) = - v ]T u n (r) v n (r)* [1 - 2 /(e„)] (3) 

n 

is the gap function that has to be self-consistently 
determined 6 , where f(e n ) = [exp(/3e„) + is the 
Fermi function with inverse temperature (3 (the sum in 
Eq. (3) being limited to positive eigenvalues e n only). 
The negative constant vq in Eq. (3) originates from the 
attractive interaction acting between fermions with op- 
posite spins (or fermionic atoms with two different in- 
ternal states) and taken of the contact-potential form 
voS(r — r'), which can be conveniently regularized in 
terms of the scattering length of of the associated two- 
body problem 9,10 . Correspondingly, the gap function has 
s-wave component only with spinlcss structure. 

Solution of the BdG equations (1) is equivalent to con- 
sidering the associated Green's function equation (in ma- 
trix form): 
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where u s = (2s+l)n/(3 (s integer) is a fermionic Matsub- 
ara frequency, 1 is the unit dyadic, and Q is the single- 
particle Green's function in Nambu's notation 11 . Solu- 
tions of Eqs. (1) and (4) arc, in fact, related by the ex- 
pression: 
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We adopt at this point Gorkov's procedure 8 for ex- 
pressing the solution of Eq. (4) in terms of the non- 
interacting Green's function Qo that satisfies the equa- 
tion 



[iw a - H(r)} g {r,r';u> s ) = S(r - r') 



(6) 



and subject to the same external potential ^(r). We are 
thus led to consider the two coupled integral equations: 
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for the normal (Gn) and anomalous (<?2i) single-particle 
Green's functions. Equations (7) and (8), together with 
the self-consistency equation 



A*(r) =v ±J2g 21 (r,v,LJ s )e- i ^ 



(9) 



(77 — > + ), are fully equivalent to the original BdG equa- 
tions (1) and (3), and hold for any coupling. 

We pass now to specifically consider the strong- 
coupling limit of Eqs. (7)-(9), and show under what cir- 
cumstances they reduce to the GP equation for spinlcss 
composite bosons with mass 2m and subject to the po- 
tential 2V(v). In this limit, the fermionic chemical po- 
tential approaches — £o/2, where e = (map)^ 1 is the 
binding energy of the composite boson which represents 
the largest energy scale of the problem. 

In the present context, achieving the strong-coupling 
limit implies that the following conditions on the poten- 
tial 

\V(r)\ « \fi\, a F |W(r)| « \V(r)\, a F \V 2 V{v)\ « \V(r)\ 

(10) 

hold for all relevant values of r. In the strong-coupling 
limit, ap ~ (2m|^|) -1 / 2 represents the characteristic 



length scale for the non-interacting Green's function of 
the associated homogeneous problem (with V(v) =0), as 
it can be seen from the expression 
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that holds for any coupling. With the conditions (10), it 
can be readily verified that the position 

&(r,r';w 8 ) ~ ^o(r - r';^| M - (V(r) + V(v'))/2), (12) 

whereby the chemical potential \i in expression (11) is 
replaced by the local form fi — (V(r) + V(r'))/2, satisfies 
Eq. (6) in the presence of the external potential. [The 
midpoint rule, albeit superfluous for most of the following 
arguments owing to the slowness of the potential over the 
distance ap, has been adopted for later convenience in 
the derivation of the current density] The approximate 
expression (12) plays in the present context an analogous 
role to the eikonal approximation for Gorkov's problem 
in the presence of an external magnetic field 8 . 

From a physical point of view, the first of the con- 
ditions (10) implies that the external potential is much 
smaller than the binding energy £ , so that the compos- 
ite boson does not break apart by scattering against the 
potential V(r). The two additional conditions (10), on 
the other hand, imply that the external potential is suf- 
ficiently slowly varying over the characteristic size a F of 
the composite boson, so that the composite boson itself 
probes the potential as if it were a point-like particle. 

In particular, for the corresponding one-dimensional 
problem with a generic (albeit non pathological) poten- 
tial V(x) that satisfies the first two conditions (10) (a_p 
therein being replaced by a — (2m|/i|) -1 / 2 ), one can 
readily verify that our position (12), namely, 
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is fully equivalent to the solution of the (one-dimensional 
version of the) Green's function equation (6) as ob- 
tained in terms of the asymptotic (a — > + ) WKB 
approximation 12 . This identification holds in the rele- 
vant range \x — x'\ < a and provided \u) s \ <§C \jj,\ (that can 
be always satisfied for sufficiently large \n\). 

Derivation of the GP equation from the BdG equa- 
tions in the strong-coupling limit exploits the position 
(12) and proceeds by using Eq. (7) to expand pertur- 
batively Eq. (8) up to third order in A(r) (the ratio 
A(r)/|/x| provides the small parameter for this expan- 
sion, by taking advantage of the "dilutcness" condition 
of the system). Equation (9) yields eventually the fol- 
lowing integral equation for the gap function: 
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with the notation 
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and 
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x So(r3,r 2 ;-a; s )^o(r3,r;w s ) . (16) 

The position (12) implies that, in strong coupling, 
Q(r, ri) vanishes for |r — ri| > ap. Since (as a con- 
sequence of the conditions (10)) A(r) is slowly varying 
over the length scale of, we can set A(ri)* ~ A(r)* on 
the left-hand side of Eq. (14) and write 
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Here, fi(r) = //— V(r) and we have introduced the bosonic 
chemical potential /j_b = 2/j+£o such that /is <C £o- Note 
that (V(r) + V(ri))/2 in the local chemical potential has 
been replaced by V(r) at the relevant order, since the 
external potential is slowly varying over the distance ap. 
Disposal of the ultraviolet divergence in Eq. (18) via the 
regularization of the contact potential 9 ' 10 in terms of ap 
and use of the strong-coupling assumption — > — oo 
have led to the approximate expressions (18) and (19). 
By a similar token, we obtain 

/"dndradra R(t, n, r 2 , r 3 ) A(n)* A(r 2 ) A(r 3 )* 
c Co |A(r)| 2 A(r)* (20) 
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Entering the approximate expressions (17)- (21) into 
the integral equation (14) for the gap function and in- 
troducing the condensate wave function 



$(r) = 
we obtain eventually: 
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(22) 



- J-V 2 $(r) + 2V(r)$(r) + ^^|$(r)| 2 $(r) = /j B $(r). 

(23) 

This is just the time-independent Gross-Pitaevskii equa- 
tion for composite bosons of mass = 2m, chemical 
potential ^b, subject to the external potential 2V{v), 
and mutually interacting via the short-range repulsive 
potential AitaE/ms (where ag = 2a p is the bosonic scat- 
tering length within the Born approximation 13,10 ). Note 
that the two-body binding energy £o has been eliminated 
from explicit consideration via the bound-state equation. 
Note also that the nontrivial rescaling in Eq. (22) be- 
tween the gap function and the condensate wave function 
has been fixed by the nonlinear term in Eq. (23). 

Equation (23) has been formally obtained from the 
original BdG equations in the limit /3/i — > — oo. In this 
respect, it would seem that this equation holds even near 
the condensation temperature T c , where the GP equation 
is instead known not to be valid. On physical grounds, 
however, at finite temperature excitations of bosons out 
of the condensate (that are not included in the present 
treatment) are expected to be important especially in 
the strong-coupling regime of the BdG equations, thus 
restricting the range of validity of Eq. (23) near zero tem- 
perature. 

The physical interpretation of the condensate wave 
function (22) can be obtained from the general expression 
for the density n(r) = (2//3) J2 S e ^P( iuj sV)Gn( r , r ; w s), 
where Gn is obtained by combining Eqs. (7) and 
(8) and expanding perturbatively up to second order 
in A(r). Recalling that the fermionic contribution 
(2/ P)^2 s ex.p(iu s r])Go{r,r;u s ) vanishes in the strong- 
coupling limit, one is left with the expression 



»(r)~-2|A(r)| : 



C/o(k,w s ) 2 £o(k, -u s ) 



^-2|A(r)| 2 (-^) =2|4>(r)| 2 



(24) 



evaluated in terms of the wave-vector representation of 
Go. Since the bosonic density ub{y) is just half the 
fermionic density n(r), from Eq. (24) we obtain that 
n B (r) = |$(r)| 2 , as usual with the GP equation. The 
normalization condition N = Jrfrn(r) (or, equivalcntly, 
its bosonic counterpart) fixes, in turn, the overall chem- 
ical potential /x (or directly ^ in the strong-coupling 
limit). 
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By a similar token, the current density can be obtained 
from the general expression 

j(r) = ^ (V - V) - J2 ^^ii(r, r'; co s )\ r=r , (25) 

' s 

by combining Eqs. (7) and (8) as before. The 
independent-particle contribution to the current now 
vanishes exactly with the midpoint rule (12). After long 
but straightforward manipulations one obtains for the re- 
maining contributions in the limit (3[i — > — oo: 

j(r) " i [m * V $W " * (r) V m * ] (26) 

which is, as expected, twice the value of the quantum- 
mechanical expression of the current for a composite bo- 
son with mass tub — 2to and wave function $(r). 

It is clear from the above derivation that higher-order 
corrections to the GP equation (23) can also be obtained 
by expanding Eqs. (7) and (8) to higher than the third 
order in A(r). In particular, to fifth order in A(r) the fol- 
lowing term adds to the right-hand side of Eq. (8) (with 
r = r') once summed over uj s : 

-|A(r)| 4 A(r)* J £ G (K ^) 3 <?o(k, -uj s ) 3 . 

(27) 

Evaluating the integral in Eq. (27) in the strong-coupling 
limit and recalling the rescaling (22), one finds that the 
term <? 3 |<I>(r)| 4 <I>(r) adds to the left-hand side of the GP 
equation (23), where the three-body interaction g% — 
— 30tr 2 a F /m is attractive. As an example, taking for mg 
the mass of 85 Rb and the typical value as ~ 250 a. u., 
one gets | <73 1 / ?i ~ 10 _27 cm 6 s _1 with the correct order of 
magnitude 14 . 

Even further higher-order corrections can be evaluated 
in this way with not much additional burden. Alterna- 
tively, when corrections to the GP equation for com- 
posite bosons are needed, one could directly proceed 
by evaluating numerically the BdG equations (1) in the 
intermediate- to strong-coupling region, thus calculating 
the physical quantities of interest in terms of the original 
fcrmionic wave functions. Such an approach would be es- 
pecially useful for addressing the intermediate (crossover) 
region, where neither the fermionic nor the bosonic char- 
acter of the system are fully developed. 

Specifically, the results of the present paper could be 
applied to determine the evolution of the density pro- 
file for a system of superfluid fermionic atoms in a trap 
when the effective fermionic attraction is increased. It 
has recently been shown 15 that, as far as the gross fea- 
tures of the density profile are concerned, this problem 
can be dealt with by solving the coupled mean-field BCS 
equations for the gap and the density within a local- 
density approximation. Such an approach reduces to the 



(bosonic) Thomas-Fermi approximation in the strong- 
coupling limit, thus missing the contribution of the ki- 
netic energy. The present derivation of the GP equation 
from the BdG equations shows, in this respect, that solu- 
tion of the BdG equations for fcrmionic atoms in a trap 
represents a refined approach that correctly treats the 
kinetic energy for all couplings. 

The connection we have established between the BdG 
and the GP equations might be also useful for deter- 
mining the limiting strong-coupling behavior of super- 
fluid fermions in the presence of an external potential, 
as, for instance, when studying strong-coupling effects in 
the Josephson and related problems. 

In conclusion, we have shown that the time- 
independent GP equation at zero temperature for com- 
posite bosons (formed as bound-fcrmion pairs) can be 
obtained on general grounds from the BdG equations de- 
scribing superfluid fermions subject to an external poten- 
tial, in the strong-coupling limit of the mutual fcrmionic 
attraction. Corrections to the GP equation have also 
been derived in powers of the condensate wave function. 

We are indebted to D. Neilson for discussions. 
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